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£-] Goostatistical A geostatistical m ethod for gradually deforming an initial 
distribution of objects, of geologic type for example, from measurements or 
observations, so as to best adapt it to imposed physical constraints o f. for example, a 
hydrodynamic type having appUcations of geostatistical modeUing of heterogeneous 
reservoirs of various objects: fracture, channels, vesicles, etc., for example. 

£-] ^_The objects bektg — are distributed in a zone of a heterogeneous medium 

according to a Poisson point process in form of figurative points with a point 
density X(x) t hat v aries according tot heir p osition ( x) i n t he z one, a r ealization o f a 
uniform random vector according to which the position of each object is defined while 
respecting density X(x) is formed, and the uniform random vector is gradually modified 
according to a gradual deformation process so as to obtain gradual migration of each 
object until a final realization best adjusted to parameters relative to the structure of the 
medium, such as hydrodynamic parameters, is obtained. 



Applicatidns : g e 




mod e lling of h e terogeneous res e rvoirs consisting of 



various obj e cts : fractur e s, chann e ls:, vesicles, etc., for exampl e . 
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BACKGROUND OF THE INVENTION 

FIELD OF THE I>JVENTIO>J 
Field of the Invention 

rOOOn ^The object of the present invention is a method for gradually deforming 

an initial distribution of objects of geologic nature, formed by simulation of a stochastic 
model of object type, from measurements or observations, so as to best adapt it to 
imposed physical constraints o f, for example, a hydrodynamic typ e for e xampl e. 

Th e method according to the inv e ntion finds applications in tho sphere of und e rground 
zonoG modelling intended to generate r e pr e s e ntations showing how a c e rtain physical 
quantity is distribut e d in a subsoil zon e (p e rm e abiUty for example) and which are best 
compatible — with obGor\^ed — or measur e d — data : — g e ologic — data; — seismic — r e cords, 
measur e ments obtained in wells, notably m e asurements of the variations with time of 
the pressure and of tho flow rat e of fluids from a reGer\^oir, e tc. 

B ACKGROU>JD OF THE R^JVENTION 

Description of the Prior Art 
r00021 French P atent FR-2,780,798 filed by the applicant assignee describes a 

method for gradually deforming a stochastic model (of Gaussian type or similar) of a 

heterogeneous medium such as an underground zone, constrained by a series of 

parameters relative to the structure of the medium. This method comprises drawing a 

number p of independent realizations (or representations) of the model or of at least part 

of the selected medium m odel from all of the possible realizations, and one or more 

iterative stages of gradual deformation of the model by carrying out one or more 

successive linear combinations of p independent initial realizations, then composite 
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realizations successively obtained possibly with new draws, etc., the coefficients of this 
combination being such that the sum of their squares is 1. 

r00031 French P atent FR-2,795,841 filed by the applicant assignee describes 

another method for gradually deforming the representations or realizations, generated 
by sequential simulation, of a non necessarily Gaussian stochastic model of a physical 
quantity z in a gridded heterogeneous medium in order to adjust them to a series of data 
relative to the structure or to the state of the medium, which are collected by prior 
measurements and observations. The method essentially comprises applying an 
algorithm allowing gradual deformation of a stochastic model to a Gaussian vector with 
N mutually independent variables, which is connected to a uniform vector with N 
mutually independent uniform variables by the Gaussian distribution function so as to 
define realizations of the uniform vector, and using these reaUzations to generate 
representations of this physical quantity z, which are adjusted in relation to the data. 

r00041 The above methods are appUcable to gridded models (pixel type models) 

suited for modelling continuous quantity fields and they are therefore ill-suited for 
modelling of zones crossed by fi-acture networks or channel systems for example, 

rOOOSI Models based on objects are spatial arrangements of a population of 

geometrically defined objects. Basically, an object type model is a Boolean model that 
can be defined as a combination of objects identical by nature with a random spatial 
distribution. Boolean models (of object type) are of great interest for the geometric 
description of heterogeneous media such as meandrous deposit systems, firacture 
networks, porous media on the grain size scale, vesicle media, etc. Geologic objects are 
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defined by their shape and size. Their location in the field is defined by taking account 
of their interactions : attraction-repulsion, clustering tendency, etc. 

fOOOSl ^Unlike pixel type models, the models based on objects can provide for 

example realistic geologic representations of an underground reservoir at an early stage 
where the data obtained by in-situ measurement are still rare. 

rOOOTI The prior art in the sphere of object type models is notably described in 

the following publications : 

-Matheron, G., 1967,— "-Element pettf — Pour une thoorie Theorie des milieux 
por e ux )) Milieux Poreux". Masson, Paris ; 

-Matheron, G., 1975,— «-"Random sets and int e gral g e om e try )) Integral Geometry". 
Wiley, New York ; 

-Serra, J., 1982, ^"Image analysis Analysis an d mathematical geology )> Mathematical 
Geology" , Vol.J, Academic Press,London ; 

-Stoyan, D.S. et al., 1995, ^"Stochastic g e ometry — Geometry and its 
applications )) A pplications", 2"** Edition, Wiley, Chichester ; 

-Lantuejoul, C, 1997, "Iterative algorithms — Algorithms for — conditional 
simulations Conditional Simulations , in Baafi and others, eds."-; Geostatistics 
Wollongong 96, Vol.1, Kluwer Acad. PubL, Dordrecht, The Netherlands, p.27-40. 

rOOOSI The position of the objects in an object type model is distributed 

according to a -the P oisson point process vvoll Icnovvn to sp e cialists . The shape and the 
size of the objects are independent of their positions. This model can be generalized by 
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a combination of objects of different nature and/or using a non-stationary density 
Poisson point process. 

r00091 Although Boolean models have been widely studied in the literature, 

there is no coherent and efficient method for constraining these models to the physical 
data, notably hydrodynamic data, which is however a major challenge for their 
application to reservoir engineering. The methods allowing gradual deformation of pixel 
type stochastic reservoir models such as those described for example in the two 
aforementioned patents carmot be directly used for Boolean models. Constraining 
Boolean models to hydrodynamic data for example requires development of coherent 
algorithms for deformation and displacement of the objects. 

SUMMARY OF THE INVENTION 

rOOIOI The method according to the invention allows te — gen e raliz e 

generalization of the gradual deformation technique described in the two 
aforementioned patents to stationary or non-stationary Boolean models, with or without 
geometric constraint as regards the well The method appears to be particularly useful 
notably for reservoir engineers anxious to coherently and efficiently adjust object type 
reservoir models. 

room The method according to the invention allows to gradually d e form 

gradual deformation of an initial realization defining the distribution of a set of objects 
in a zone of a heterogeneous medium such as a geologic structure, formed by simulation 
of an object type stochastic model, the objects being distributed in the zone according to 
a Poisson point process in form of figurative points with a point density A,(x) that varies 
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according to their position (x) in the zone. It essentially comprises the following-4we 
stages -steps : 

-generating a realization of a uniform random vector according to which the position of 
each object is defined while respecting density X(x)^; and 

-gradually modifying the imiform random vector according to a gradual deformation 
procedure, so as to obtain gradual migration of each object and consequently gradual 
change in the distribution of the objects in the zone, until a final realization best 
adjusted to parameters relative to the structure of the medium, such as hydrodynamic 
parameters, is obtained. 

f001 21 It is possible to limit the migration of the figurative points in a 

subdomain of the zone (a well crossing the zone for example) by imposing a zero point 
density in the complementary part of the subdomain. 

r00131 A ccording to an implementation mode, one gradually changes from a 

realization containing a first set of Ni points to a realization containing a second set of 
N2 points by constructing a chain N(t) of Poisson numbers between the two numbers Nj 
and N2, using the gradual deformation procedure. 

f00141 _It is also possible to gradually deform the size, the shape and the 

orientation of an object during its migration by using the gradual deformation 
procedure. 

r00151 In cases where it contains errors, it is also possible to gradually adjust the 

point density A.(x) by using the gradual deformation procedure. 
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r00161 The final realization obtained by means of the method is a realistic 

representation of the configuration of the objects in the heterogeneous medium. 

r00171 The method according to the invention finds applications in the sphere of 

underground zones modelling intended to generate representations showing how a 
certain physical quantity is distributed in a subsoil zone (permeability for example) and 
which are best compatible with observed or measured data : geologic data, seismic 
records, measurements obtained in wells, notably measurements of the variations with 
time of the pressure and of the flow rate of fluids firom a reservoir, etc. 



BRffiF DESCRIPTION OF THE DRAWINGS 

r00181 Other features and advantages of the method according to the invention 

will be clear firom reading the detailed description hereafter, with reference to the 
accompanying drawings wherein : 

r00191 [-1 Figur e s Figs. 1 A. IB respectively show a density function of a non- 
stationary Poisson point process and a realization of the non-stationary Poisson point 
process generated by the sequential method ; 

[00201 [-1 Figures F igs. 2A to 2D show various examples of gradual migration 

trajectories between two points of a non-stationary Poisson point process ; 

r002n [-1 Figures Figs. 3A to 3H show various successive stages of a 

realization of the non-stationary Poisson point process dviring the gradual migration of 
the points ; 
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r00221 [-] Figur e Fig. 4 shows the possible domains of migration of a disk in 

the case of three conditioning points ; 

r00231 [-1 Figur e Fig. 5 shows, by way of example, a complete chain of 

successive realizations of a Boolean simulation of elliptical objects where deformation 
parameter t ranges between -tt and n by intervals At=0.l7i ; and 

r00241 [-1 Figur e F ig. 6 shows, by way of example, an incomplete chain of 

successive realizations of a Boolean simulation of elliptical objects where deformation 
parameter t ranges between 0 and 0.2n by intervals At=0.0l7r. 

DETAILED DESCRIPTION 

General points 

r00251 ^The geologic objects to which the method applies are for example more 

or less large fractures inside a reservoir zone, or channels. The method can also apply to 
granular or vesicle type structures of much smaller size. All these objects are difficult to 
model by means of pixel type models. 

r00261 The progressive migration operations described hereafter apply to an 

initial model where the positions of the objects are represented by point configurations 
(referred to as point procedures). The distribution of these points varies in density 
according to their positions in the modelled zone. This distribution is based on various 
data known from measurements or observations : geomechanical measurements 
obtained in wells for example, seismic data obtained from prior seismic operations. 
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[00271 Starting from this initial model, the method allows to gradually doform 

gradual deformation of the initial distribution according to a series of rules so that, in 
the final distribution, the model is optimized in order to better correspond to constraints 
imposed on one or more physical parameters such as, for example, a distribution of 
permeability values. The migration rules are such that a great number of different points 
of the model can be globally displaced from a limited number of control parameters. 

r00281 An objective function measuring the difference between the physical 

parameters from the real heterogeneous medium and those simulated on a realization of 
the stochastic model is generally defined. The value of the objective function therefore 
depends on control parameters of the stochastic model. The values of these control 
parameters are obtained by minimizing the objective function. 

Summary of the Poisson point process 

r00291 The Poisson point process is a coimtable random set of points distributed 

in any space This set of points has the following characteristics : 

f00301 Let D be a domain of 5H". If the volume of D, denoted by I D I, is finite, 

then the numbers of points lying in D, denoted by N(D), follows a Poisson' s law of 
parameter X Id I. Let : 

p[NiD) = nY'^'^ ^^J^ V« > 0 (1) 

where X is called the point process density ; it measures the average nvimber of points 
lying in a unitary volume domain of 5R". 
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rOOSn Let Di, D2,..., Dk be domains of 9?° unconnected two by two, then the 

numbers of points lying in these domains N(D,), N(D2),..., N(Dk) are mutually 
independent random variables. 

[00321 Conditionally on N(D)=np, these np points are independent and evenly 

distributed in D. 

r00331 The description hereafter presents the Poisson point process in the finite 

domain D. 

Migration of a stationary Poisson point process 

r00341 Consider the problem of the migration of a realization of a stationary 

Poisson point process in D rectangular. In order to simplify the presentation, D is 
assumed to be a unitary hypercube [0,1]" with n dimensions. Let Xi and X2 be two 
independent points uniformly drawn from [0,1]". We define a trajectory between and 
X2 by : 

x{t) - g[g-' (jc, ) cos ^ + (x^ ) sin r J (2) 

where G is the centered and reduced Gaussian distribution function. According to the 
gradual deformation algorithm described in the aforementioned French p atent ¥Rr- 
2,780,798, for any t, x(t) is a uniform point in [0,1]". When the two points are fixed, the 
trajectory of the gradual migration between them is entirely determined. The place 
change of one of the two points will change the migration trajectory. We-It can be 
shown that the trajectory defined by Equation (2) is symmetrical in relation to the centre 
of the domain [0,1]", whatever the number of dimensions n. This suggests that, even if 
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the two points are isolated in a comer/side of the domain, the trajectory of the gradual 
migration between them can always reach the opposite part of the domain. 

Migration of a non-stationary Poisson point process 

r00351 Consider now the Poisson point process in a domain D of density X(x) of 

general shape. The number of points in D is a random Poisson variable of average 
A.(D) = j^A(x)dx . These points are independently distributed in D according to the 

probabiUty density function : 

f(x) = Mx) I ?i{P) x^D (3) 

f00361 Simulation of a Poisson point process of density in D can be carried 

out in two stages : 

- generating a number n according to Poisson' s law of average equal to X(D), then 

- generating n points in D irrespective of one another according to the same probability 
density f(x). 

[00371 The example of Fig. 1 A shows a density function, and Fig. IB shows a 

realization of a Poisson point process of density ?i(x). 

r00381 If law f(x) is simulated by inversion of its distribution function, then 

point X corresponds to a uniform vector u. The gradual migration algorithm can thus be 
applied to the Poisson point process of density A,(x). By construction, this method 
preserves the density and the number of points of the initial process. Figs. 2A to 2D 
show four migration trajectory examples. 
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[00391 Take the example of the non-stationary Poisson point process in a two- 
dimensional domain D = [0,1]^ whose density increases linearly in axis X and is 
constant in axis Y. Let (x,y) be the vector of the coordinates of a point in D of this 
process. Then, (x,y) admits the bivariable probability density : 

f{x,y) = 2x, (^,j)e[0,lf (4) 

r00401 Simulation of a point according to the above law (4) is simple : 

- generating abscissa x according to the linear density f(x) = 2x, then 

- generating ordinate y uniformly between 0 and 1 . 

r004n Let : 

X = Vw 

y = v (5) 

where u and v are two independent and uniform numbers between 0 and 1 . 

r00421 W e can thus apply Thus a g radual migration algorithm (2) can be applied 

to the uniform vector (u,v) in order to establish a migration trajectory for point (x,y) in 
D. 

r00431 In practice, density function ?i(x) often comes in form of a grid. By way 

of example, consider the case of a two-dimensional grid of M x N nodes. Let Xj and yj 
be the coordinates of node (i j). The marginal law for Xj is then : 

f{x,) = j^f{x,,yj) (6) 
and the conditional law for yj, knowing Xi, is : 

f.^yj) = Kx,,yj)ifix,) (7) 
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r00441 The points of the Poisson process of density X(x) can thus be readily 

generated on a grid from these monovariable probabilty laws. 

r00451 ^By w ay o f i Uustration, we have constructed a t wo-dimensional d ensity 

function is constructed from a centr e d centered and reduced Gaussian simulation. The 
variogram is anisotropic and of Gaussian type. The principal anisotropy directions are 
diagonal in relation to the coordinates system. The scale factors in the principal 
anisotropy directions are 0.3 and 0.1 respectively. The size of the field is 1x1 and it is 
discretized in 1000x1000 pixels. The Gaussian numbers are changed into positive 
numbers according to the following expression : 

= 400^^^'^ (8) 

where Y(x) is the Gaussian simulation. Fig.lA shows the density function thus 
constructed. Fig. IB shows a Poisson point process realization admitting the density 
function of Fig.lA. 2000 points are generated according to the sequential method. 

r00461 ^Figs. 3A to 3H illustrate the evolution of a Poisson process realization 

during gradual migration of the points. It can be seen that the point density is respected 
during the migration. 

Migration of conditional Boolean simulation objects 

[00471 ^The above method can be immediately applied to the migration of a 

domain S of any form in D. In fact, the migration of a point in S can be carried out in D 
using the truncated probability density function : 

f,ix) = Mx)l^^, I ?i{S\ x^D, S^D (9) 
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[00481 ^Using the Markovian iteration algorithm described in the reference to 

Lantuejoul, 1997, it is possible to simulate a Boolean model in a domain D knowing 
that two subsets Ci and Co of D respectively belong to the union of the objects and to its 
complement. Then, the migration algorithms in any domain can be used for gradual 
deformation of the Boolean simulations conditioned by geometric well data. In fact, 
from a conditional Boolean simulation and without compromising conditioning by Ci 
and Co, the objects must move only in their respective domains defined according to the 
shape of the objects and the configuration of Ci and Cq. 

[00491 Consider an object A of a conditional Boolean realization, which 

includes a subset Cia of Ci and excludes C©. If, during migration of A, it still has to 
include Cia but excludes Co, the allowed migration domain for object A is : 

D, = {x:C,nA^=0;C,,czA^} (10) 

[00501 ^Fig -4 shows a case with three conditioning points and the eight possible 

domains of migration of a disk. If, for example, the disk is supposed to always cover 
points (a) and (b) but avoid point (c), then its centre can move only in domain 3. 

Migration with appearance and disappearance of points 

[00511 The number of points in D of a point process follows a Poisson's law of 

parameter X(D), It is therefore necessary to vary the number of points in D during their 
migration. In this section, we first pres e nt two methods of constructing Poisson number 
chains are first presented , and then w e shall s ee h ow to migrate between two sets of 
points whose cardinals are not identical is presented . 
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r00521 We tr\^ to construct a A chain of Poisson numbers is attempted to be 

constructed b etween two Poisson numbers Ni and N2, independently generated by 
inversion of the distribution function. Let U| and U2 be two independent uniform 
numbers (between 0 and 1) from which numbers Ni and N2 are obtained. 

N,^F-\U,) 

N,=F-\U,) (11) 

where F ' represents the inverse distribution function of Poisson' s law. According to the 
gradual d eformation a Igorithm, we con construct a c hain o f u niform n umbers can be 
constructed b etween Ui and U2 by : 

U(t) = G[G-'(U,)cost + G-\U^)smt\ (12) 

r00531 ^By inversion of the distribution function of Poisson' s law, we obtain a 

chain of Poisson numbers is obtained . 

Nit) = F-'[U{t)] (13) 

r00541 The inverse distribution function can be calculated by dichotomy. 

However, if the parameter of Poisson' s law is too great, this method remains quite 
costly. 

r00551 In order to save calculating the inverse distribution function of Poisson' s 

law, another way of generating Poisson's numbers can be e nvisaged considered . It is 
well-known that the number of jumps of a Poisson process of parameter 1, in a segment 
of length X, precisely follows a Poisson's law of parameter X, Simulation of a Poisson 
process is easy to do. In fact, the intervals between two consecutive jumps of the 
process are independent and follow the exponential law of parameter 1. By gradually 
deforming the exponential segments, wo obtain a chain of Poisson numbers is obtained . 
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It is easy to construct a chain of exponential numbers because calculation of the inverse 
distribution function of the exponential law is easy. It can be reminded here that the 
distribution function of the exponential law of parameter 1 is written as follows : 

F{s) = l-e-\ s>0 (14) 
and its inverse function is : 

F'* (m) = - ln(l - w), 0 < w < 1 (15) 

r00561 The fact remains that, since the number of jumps in the segment of length 

X. is X on average, the number of chains of exponential numbers is close to X on average. 

r00571 Consider now N ow the migration of a set of Ni points (set 1) to a set of 

N2 points (set 2) is considered . Since Ni ^ N2, the migration from set 1 to set 2 
necessarily implies the appearance or the disappearance of certain points. The number 
of points that will appear or disappear is determined by the chain of Poisson numbers 
between Ni and N2. The migration algorithm is as follows : 

a) Calculate the maximum number N^ax of chain N(t). 

b) Complete set 1 with N^ax - Ni points and set 2 with N^x - N2 points. 

c) For each set, order the points from 1 to Nmax- AH the initial points are ordered 
from the beginning. 

d) Calculate the migration trajectory of point n of set 1 to point n of set 2 
(n=l,2,...,N^x). 

e) For each state t of the set, remove the N^ax - N(t) last points. 
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CLAIMS 

1. (Currently Amended) A method for gradually deforming an initial 
realization formed from measurements or observations and defining tke-a.distribution of 
a set of objects in a zone of a heterogeneous mediu m such as a g e ologic otructure , 
generated by simulation of an object type stochastic model, the objects being distributed 
in the zone according to a Poisson point process in a_form of figurative points with sl 
point density X.(x) varying according to their position (x) in the zone, characterized in 
that it compris e s comprising : 

I'l generating a realization of a uniform random vector according to which fee-a 
position of each object is defined while respecting density X(x)^i and 
£-1 gradually modifying the uniform random vector according to a gradual deformation 
procedure, so as to obtain gradual migration of each object and consequently gradual 
change in the distribution of the objects in the zone, until a final realization best 
adjusted to parameters relative to the structure of the medium, such as hydrodynamic 
paramotors, is obtained, which gives a realistic representation of the configuration of the 
objects in the modelled heterogeneous medium. 

2. (Currently Amended) A method as claimed in claim 1, charact e rized in that 
wherein m igration of the figurative points representing objects in a subdomain of the 
zone is limited by imposing a zero point density in the complementary part of the 
subdomain. 

3. (Currently Amended) A method as claimed in claim 1, characterized in that 
ono gradually changoG from wherein a realization containing a first set of Nl points is 
gradually changed t o a realization containing a second set of N2 points by constructing 
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a chain N(t) of Poisson numbers between the two numbers Nl and N2 using the gradual 
deformation procedure. 

4. (Currently Amended) A method as claimed in claim 1, charact e rized in that 
tbe -wherein size, Sie-shape and the-orientation of an object are gradually modified 
during ite -it's the migration of the object using the gradual deformation procedure. 

5. (Currently Amended) A method as claimed in claim 1, charact e riz e d in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

6. (Currently Amended) A method as claimed in claim 2, charact e riz e d in that 
one gradually changes from wherein a realization containing a first set of Nl points is 
gradually changed to a realization containing a second set of N2 points by constructing 
a chain N(t) of Poisson numbers between the two numbers Nl and N2 using the gradual 
deformation procedure. 

7. (Currently Amended) A method as claimed in claim 2, characterized in that 
the -where size, the-shape and the-orientation of an object are gradually modified during 
its -it's the m igration of the object u sing the gradual deformation procedure. 

8. (Currently Amended) A method as claimed in claim 3, characteriz e d in that 
the -wherein size, tfee-shape and Sie-orientation of an object are gradually modified 
during its -it's the m igration of the object using the gradual deformation procedure. 
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9. (Currently Amended) A method as claimed in claim 2, charactorizcd in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

10. (Currently Amended) A method as claimed in claim 3, characterized in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

11. (Currently Amended) A method as claimed in claim 4, charact e rized in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

12. (Currently Amended) A method as claimed in claim 6, characterized in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

13. (Currently Amended) A method as claimed in claim 7, characterized in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure. 

14. (Currently Amended) A method as claimed in claim 8, charact e rized in that 
wherein p oint density ((x) is gradually adjusted using the gradual deformation 
procedure 
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15. (New) A method as claimed in claim 1 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

16. (New) A method as claimed in claim 2 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

17. (New) A method as claimed in claim 3 wherein the mediiun is a geologic 
structure and the parameters are hydrodynamic parameters. 

18. (New) A method as claimed in claim 4 wherein the mediiun is a geologic 
structure and the parameters are hydrodynamic parameters. 

19. (New) A method as claimed in claim 5 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters, 

20. (New) A method as claimed in claim 6 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

21. (New) A method as claimed in claim 7 wherein the medivmi is a geologic 
structure and the parameters are hydrodynamic parameters. 

22. (New) A method as claimed in claim 8 wherein the medium is a geologic 
structxire and the parameters are hydrodynamic parameters. 
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23. (New) A method as claimed in claim 9 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

24. (New) A method as claimed in claim 10 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

25. (New) A method as claimed in claim 11 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

26. (New) A method as claimed in claim 12 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

27. (New) A method as claimed in claim 13 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 

28. (New) A method as claimed in claim 14 wherein the medium is a geologic 
structure and the parameters are hydrodynamic parameters. 
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